The linear spectrum of a quantum dot coupled to a nano-cavity 
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We develop a theoretical formalism to model the linear spectrum of a quantum dot embedded in 
a high quality cavity, in presence of an arbitrary mechanism modifying the homogeneous spectrum 
of the quantum dot. Within the simple assumption of lorentzian broadening, we show how the 
known predictions of cavity quantum electrodynamics are recovered. We then apply our model to 
the case where the quantum dot interacts with an acoustic-phonon reservoir, producing phonon 
sidebands in the response of the bare dot. In this case, we show that the sidebands can sustain 
the spectral response of the cavity-like peak even at moderate dot-cavity detuning, thus supporting 
recent experimental findings. 
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I. INTRODUCTION 

The description of semiconductor quantum dots (QDs) 
as isolated atomic-like quantum systems is largely over- 
simplified. The confined electrons and holes in a QD, in 
fact, interact rather efficiently with both the electronic 
and vibrational [ll-IH degrees of freedom of the semicon- 
ductor environment, in ways that can only be described 
beyond the simple perturbation theory. In addition to 
the semiconductor medium, QDs interact with the sur- 
rounding electromagnetic field, especially if embedded 
in a photonic structure with sharp electromagnetic reso- 
nances. As an example, the electromagnetic field can vc- 
hiculate an excitation transfer between two distant QDs 
with non-overlapping electronic states [B-Q ■ If the QD is 
embedded in a high-quality cavity, the 3D confinement of 
electromagnetic field can lead to observation of the strong 
coupling between one QD and the resonant mode of the 
electromagnetic field |9l4l3l|. This system, however, can 
not be seen as a perfect parallel to the atom-cavity cou- 
pling in cavity quantum electrodynamics (CQED) 
The semiconductor environment, in particular, can af- 
fect the system in several ways that have no analogous 
in its atom-cavity counterpart. A first effect is brought 
by the coupling of the QD to an external reservoir (e.g. 
phonons) that can produce a significant change in the ho- 
mogeneous spectral signature of the dot. As examples we 
quote the broad sidebands originating by the coupling to 
longitudinal acoustic phonons beyond perturbation the- 
0, H, [sl. [isl- 171. or the similar effect due to opti- 
19|. In addition to these homogeneous 
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cal phonons 

modifications of the bare QD spectrum, other significant 
spectral changes can arise when the QD is multiply ex- 
cited, due to transitions between continuum states in the 
wetting layer above the QD confining barrier. This effect 
occurs already at moderate excitation and leads to a size- 
able enhancement of off-resonance light emission j20l - |22j . 

Here, we address the first class of these semiconductor- 
related effects, where a homogeneous change in the spec- 
trum of the single QD is present. We restrict to the linear 
spectral response of the cavity-QD system, holding at low 



excitation density. Using a Maxwell formalism, we show 
how this effect translates in the overall spectral signature 
of the QD-cavity system, at varying detuning and QD pa- 
rameters. Our model accounts for the specific shape of 
the resonant cavity mode as well as for the microscopic 
parameters of the QD excitonic transition. 

We first develop the general formalism, assuming an 
arbitrary energy- dependent self-energy for the QD exci- 
tonic transition. Then we apply the model to two cases. 
The first is that of a lorentz-shaped QD line, character- 
ized by a constant broadening 79. In this case, we recover 
the result that is well known from CQED [2^ . Then, we 
assume a QD spectrum arising from the coupling to lon- 
gitudinal acoustic phonons, with the corresponding self- 
energy modeled within a second-Born approximation 0] . 
The coupling produces phonon sidebands in the bare QD 
spectral response. We show how the spectral signature of 
these sidebands is enhanced by the presence of the cavity 
resonance, even for moderate QD-cavity detuning. The 
persistence of light emission at the cavity-like peak has 
been the object of several experimental investigations re- 
cently [1^ |20 22,, ■ small detuning, the acoustic 
phonon mechanism is expected to contribute significantly 
to this effect, as the energy width of the sidebands is de- 
termined by the exciton spatial confinement [l5| . and 
amounts to 1-2 meV in typical samples. This effect has 
been recently addressed using phenomenological dephas- 
ing to account for the modified QD spectral signature 
or within a more microscopic approach to the 



phonon sideband mechanism at zero QD-cavity detuning 
jl7| . The importance of our work lies in the fact that a 
general homogeneous mechanism acting on the QD spec- 
trum is modeled, and an explicit expression for the total 
emission spectrum, accounting for spatial and spectral 
cavity form factors, is derived. 

In Section II, we present the general theoretical formal- 
ism. Section III is devoted to deriving the simple CQED 
result in the limit of lorentz QD broadening. In Section 
IV, we study the case of a QD coupled to a reservoir of 
longitudinal acoustic phonons, and discuss how phonon 
sidebands enhance the cavity-like emission spectrum at 
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finite cavity QD detuning. In Seetion V, we present our 
conclusions. 



II. THEORETICAL FORMALISM 

A. Maxwell equations 

We consider the system of one QD embedded in a reso- 
nant nanocavity. The cavity can be of any kind (e.g. pil- 
lar photonic crystal defect [T]|. microdiscs[l3|. etc.), 
with the only assumption that one well distinct resonant 
mode exists in the vicinity of the QD transition wave- 
length. Our objective is to derive the physical parame- 
ters characterizing an effective coupling to this mode, at 
frequency Uc, from the microscopic details of the electro- 
magnetic field in the structure. We assume a QD lying 
at position ro. Typically, this position is selected to lie 
where the electric field has maximum amplitude, in order 
to maximize QD-cavity coupling. In the limit of low QD 
excitation, the spectra are determined by the linear op- 
tical response, and are described by Maxwell equations 
for the electric field coupled to the linear susceptibil- 
ity tensor of the QD. Under this assumption, the steps 
leading to a set of coupled mode equations are formally 
the same as in our previous works 0, 0] ■ In particular 
Maxwell equations are cast into an integral Dyson equa- 
tion [2^. We denote with e = e(r) the spatially depen- 
dent dielectric constant that characterizes the resonant 
photonic structure. In the frequency domain, we have 
(assuming non magnetic medium and no free charges) : 



V A V A£(r,w) 



,2 r 



e{r)£ (r,w) 



+47r / dr'xQD {r, r',uj) -S {r',uj) 



where £ {r, lu) is the electric field, r is the 3-D position 
vector, and xqd the 3x3 linear optical susceptibility 
tensor of the QD subsystem. In order to define a hermitic 
problem, we adopt the standard replacement [29j 



and 



T = 




This leads to the following hermitic problem: 



(1) 



(2) 



TQ(r,w)-^Q(r,c^) 




dr'xQD {r, r',uj) 



(3) 



B. Photon Green's function 

We introduce the in-plane Green's tensor of the pho- 
ton 0{r,r',uj), which is defined as the Green's tensor of 
the Maxwell equation. We have previously shown that 
simple analytical expressions hold in the case of a QD in 
a homogeneous medium @ or in a planar microcavity ■ 
In the general case, a compact analytical expression can- 
not be found. Formally, the Green's function is defined 
as : 



U! 



T(r) 



Q{r, r' ,Lj) = 5{r — r') 



(4) 



where T(r) is a time independent, hermitian, linear dif- 
ferential operator that possesses a complete set of eigen- 
functions {$„(r)} where u is a continuous index. The 
set is considered as orthonormal. This differential prob- 
lem belongs to the class described by Fredholm theory 
[sot ■ It therefore admits a formal solution in terms of the 
resolvent representation 



du 



(5) 



Once obtained the Green's function of the photonic struc- 
ture, the solution of Eq. ([3]), corresponding to an input 
field Qo('", w) can be written as follows j28j : 



(6) 



-47r 



— dr dr g[r,r ,uj)—^ — == Q{r,uj), 

^ J J 



The key assumption of our procedure is that one strongly 
resonant mode exists and is energetically well distinct 
from any other spectral feature (discrete or continuous) 
of the structure under investigation. This is the case for 
all kinds of high-quality nanocavities. Close to resonance 
LU M ujc, the following approximation then holds 



g{r, r',uj) 



2iOc{uJc — cj — i^) 



du 



2a;„(a;„ - uj - i^) 



(7) 



A similar expression was used by Sakoda et al. [3l| and 
Hughes et al. [s^] ■ Here, we neglect the longitudinal 
optical modes, consistently with the exciton optical se- 
lection rules that we assume (see below). In compact 
form we obtain 

g(../,.). *°'-'^s"-'":' (8) 



if) 



The resonant cavity mode arises as sharp resonance in the 
energy-dependent density of the eigenmodes. We have 
characterized this resonance by a damping constant k. 
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that models the finite hfetime of the mode. This step is 
necessary, as we are approximating an everywhere con- 
tinuous mode spectrum with one discrete mode plus a 
nonresonant continuum. Formally, this passage can be 
justified in terms of the quasi-mode theory p3L Is^ ]. by 
assuming weak coupling between an ideal undamped cav- 
ity mode and the vacuum electromagnetic field outside 
the cavity. In Eq. (|5]), gc{r,r') represents the contri- 
bution of all other modes, and is supposed to be small 
at w « Wc- A complete numerical calculation of cavity 
eigenmodes, like e.g. that carried out in Ref. (fssh. can 
be used to test this assumption. In the following, we will 
express gdi^, r') as the sum of its real and imaginary parts 
a{r,r') and ib{r,r'). As shown later, these are responsi- 
ble - respectively - of a shift and a broadening of the 
QD emission spectrum. More precisely, the term 5(r, r') 
is responsible of the decay of the excited QD into the 
continuum of background electromagnetic modes. This 
determines the free decay rate of the QD, usually denoted 
as 7 in CQED. 



III. CQED LIMIT 

As a simple test of our formalism, we can recover the 
limit of one two-level emitter in a resonant cavity, namely 
the simplest CQED system. Our derivation has the ad- 
vantage of relating all CQED parameters to microscopic 
expressions for the semiconductor QD - nanocavity sys- 
tem under investigation. 



A. QD susceptibility tensor 

In semiconductors with cubic symmetry (e.g. InGaAs), 
the QD susceptibility tensor is expressed as 



XQD (r, r', c^) = ^""^ ^""^ 



1 
1 




(9) 

where Hcv is the Bloch part of the interband dipole matrix 
element and 5'(r) is the electron- hole wave function in the 
QD, taken at r = re = r/j. Here, we are assuming heavy 
holes only, hence the z-component is uncoupled to the 
electromagnetic field. We further assume to deal with 
a single QD transition, having one specific polarization 
(e.g. along x). Then, the susceptibility tensor is replaced 
by a scalar, where 



XQD (w) = 



1 



LUq 



(10) 



Here, 70 is an additional non-radiative damping rate of 
the bare QD resonance. Given the small size of the QD 
with respect to the cavity mode spatial extension, we can 
safely approximate 4' (r) = (5(r — ro). Then, Eqs. ^ and 



(HI) result m 

Q{r,uj) = Qo (r,w) + 



M (w) 



2wc(wc - w - if) 



with 



M{uj) 



hc^ (wq — oj — i^) 



(11) 
(12) 



where em = e(''o). This expression is the starting point for 
computing the spectral properties of the cavity-QD system. 
It gives direct access to the linear response spectrum of the 
system. From this, the emission spectrum can also be mod- 
eled. 



B. Emission spectrum 

We first take all fields at position ro. Then, Eq. can 
be rewritten in compact form 



Q (r(,,aj) ((oJo - co- i^)(cjc - u) ~ i'^) - g^) 
= Qo {ro,i^) (a;o — a; — i^)(cjc — oj — i^) , 



(13) 



with 



7 = 70 + 



It = 



h€M 



(14) 
(15) 
(16) 

(17) 



We then compute the emission spectrum of the system from 
the linear response equation (|lip . using the virtual oscillating 
dipole method [s^. The method is based on the assumption 
that spontaneous emission is the linear response of the system 
to vacuum field fluctuations. We therefore solve Eq. (|lip with 
Qo (r) given by the field produced by an oscillating dipole at 
the QD position, in the photonic structure. We define: 



Sq (oj) 



{ujo —uj — i^){idc — uj — if) — g|2 
Using equation ([ISp . we find: 



(18) 



Q (r, oj) = Qo (r, uj) + 



47rcj ^i^„ 



Sg (uj) 



*o(r)$S(n))c^ 



2llIc{u}c — (jJ — : 



+ gc{ro,ro) 



Qo {ro, 



(19) 



The input field Qo (r, co) in our formalism is the field present 
in the photonic structure in the absence of the QD. This field 
can be computed by a Green's function procedure similar to 
the onepresented above. In this case, following the method 
of Ref [2g|, the background dielectric system is the free space. 
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while the perturbation is the photonic structure itself. This 
procedure is presented in Appendix A. We further neglect the 
first and third term on the right-hand side of Eq. 1191 as they 
are off-resonant with respect to the cavity mode. We obtain 
the following expression for the emitted field 



Ve{r)S {r, oj) = 
Anuj^lJ,^^ Sq (lj) 



$o(r)$S(''o)c^ 



So (royi^) 



(20) 

This result can be now traced back to the well known ex- 
pressions for the QD-cavity emission spectrum [stI . |38| . We 
use the fact that Qo (r, uj) is in general smoothly varying as a 
function of uj, as discussed in Appendix A. Hence, we assume 
£o {ro,u)) £o. Then 



So 



(21) 



with the semiconductor cavity form factor expressed as 



(7-«)^ 



'l>o(r-)*S('^o)c^ 



(22) 

The remaining factor in Eq. (|21[) is the emission spectrum of 
the QD, as in atomic CQED, 5 (cj) = \Sq (oj) p, expressed in 
the resonant case {ujc = i^o) as 



5(^) = 



UJ — + 



(23) 



with 



-ujo - 7(7 + K)d 



(24) 



This is the usual CQED result [13, H^. We see from equa- 
tion p4|) that the coupling of the QD to the electromagnetic 
field of the modes other than the cavity mode, produces a 
radiative shift and an additional radiative damping, respec- 
tively proportional to the real and imaginary parts a{ro,ro) 
and b{ro, ro) of the photon Green's function. The shift simply 
redefines the resonant frequency and will be neglected in the 
following. The background electromagnetic field however, has 
also an impact on the radiative damping of the QD. In partic- 



ular, 7r 



originating from the term b{ro,ro) 



models the radiative decay into the nonresonant background 
electromagnetic field. Within the CQED formalism, the QD 
decay rate denoted by 7 can be linked to our result by defining 
the damping rate 7 = 70 4- 7r ■ 

For example, we use Eq. (|14|) to compute the Rabi split- 
ting of a semiconductor QD embedded in a photonic crystal 
nanocavity. We model the optical cavity mode as a Gauss- 
shaped mode ^oif) with spatial extension corresponding to 
the typical size of a mode in this system. We further assume 
zero QD-cavity detuning, and the QD position centered at 
ro = 0. By defining the volume Vm of the Gauss mode, we 
find: 

/ . . (25) 

This expression coincides with that obtained in Ref. [stI . With 
realistic numerical values for InAs QDs in photonic crystal 
nanocavities (/i^„ = 4S>QmeVnm?, Vm = 0.04/im"^), we find 
g « 200^16 V". 



IV. BEYOND THE MACROATOM PICTURE 

Recent studies have demonstrated that a semiconductor 
QD displays spectral features that depart from the simple 
picture of a two-level system. One typical example is the 
non-perturbative coupling to acoustic phonons, resulting in 
broad phonon sidebands in the exciton spectrum. This mech- 
anism has now been extensively characterized both theoret- 
ically [3, i, m and experimentally 0, Q, 01 ■ Another mech- 
anism that has been recently investigated is the transition 
between multi-exciton manifolds, involving the continuum of 
excited states of each manifold (sometimes referred to as 
"shakeup process"). This mechanism has proven very effec- 
tive especially when a QD is embedded in a resonant cavity, 
giving rise to intense PL at the cavity mode even at very large 
cavity-QD detuning - the cavity feeding mechanism. The for- 
malism discussed here can be generalized to situations like 
the first one, characterized by a homogeneous spectral modi- 
fication, by replacing the simple QD susceptibility (|10|l with 
the appropriate model. Here, as an example, we discuss the 
case of exciton-acoustic phonon coupling with formation of 
phonon sidebands. We are still interested in determining the 
emission spectrum in the form (|21[l. 



A. QD susceptibility tensor 

The coupling of one exciton to the LA-phonon band is de- 
scribed exactly, through the solution of the independent Boson 
model [l^,|43l. It has however been shown 0] that a very good 
account of the exciton spectrum can be obtained already at 
the 2nd Born perturbation level, with the advantage of hav- 
ing a simple expression for the exciton-phonon self-energy. 
We thus rewrite the QD exciton susceptibility including the 
exciton-phonon self energy as 



XQD {uj) = 



ujo-uj-i^ + T.[uj) 



(26) 



where, within second Born approximation and restricting to 
only one phonon band. 



s(^) = E 



-uj{q) 



9% \\n(si)) 



UJ ^i- 



- UJO ■ 



UJ + i- 



' uJo 4- uj{q) 



Here, n(g) is the Bose-Einstein equilibrium phonon occupa- 
tion at temperature ksT . We consider the case of deforma- 
tion potential coupling with acoustic phonons of dispersion 
ujn = gs (? =1 f? I), where s is the sound velocity, as in Ref. 
[15I . In Fig. [T]we display the imaginary part of the QD sus- 
ceptibility, as computed at ksT = \QK for an InAs QD of 
lOnm diameter. It should be noted that the phonon spectral 
features do not depend specifically on the shape of the exciton 
wave function but are only determined by its volume [l^. In 
the plot, we notice the pronounced sidebands compared to the 
spectrum of an ideal exciton. The sidebands are more pro- 
nounced on the high energy side, where they are determined 
by acoustic phonon emission. 



B. Emission spectrum 

Intuitively, the emission intensity at the cavity-mode fre- 
quency depends on the optical density of the underlying ex- 
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FIG. 1: Imaginary part of the quantum dot susceptibility in 
presence of LA-phonon coupling (full) and without phonons 
(dashed), computed at T=10K. As an illustration, we plot 
the cavity mode optical density at 1 meV positive detuning 
(dotted). 



citon spectrum. Hence, the presence of acoustic phonon side- 
bands is expected to enhance this PL intensity, when the cav- 
ity is detuned from the exciton. This is illustrated in Fig. [1] 
where the cavity mode spectrum is plotted at 1 meV positive 
detuning with respect to the exciton peak. 

We use the QD susceptibility 1)26^ to compute the emission 
spectrum (|21|l . The form factor J^{r,oj) is still expressed as 
while the QD emission spectrum now reads 



(7-«)" 



{u!o — U! — + E(tj))(cjc — LJ — i^) — g'^ 



(27) 



As expected, the exciton-phonon coupling results in a mod- 
ified emission spectrum. In particular, the exciton-phonon 
self energy is responsible for a modified intensity at the cav- 
ity mode frequency and a small polaron shift of the exciton 
frequency. In Fig. (2] we plot the computed spectrum at 
fcflT = 10 K for various values of the exciton-cavity detun- 
ing. While the strong coupling features remain essentially 
unchanged for zero-detuning (see panel b), we can clearly see 
in panels (a), (c), and (d) that phonon sidebands can effi- 
ciently emit through the cavity mode, as also found by other 
theoretical approaches [2^. l4l|. As a consequence, the peak 
at frequency « cjc is enhanced with respect to the simple 
CQED model, provided the detuning is not larger than the 
energy extent of the phonon bands. This effect was widely 
investigated experimentally during the last years. It has been 
reported for QDs embedded in various systems such as pho- 
tonic crystal nanocavities [l^, [2ll . and micropillars [24|. 
We can see in Fig. [2l^d) that for detuning exceeding the typ- 
ical broadband width, this feature starts to disappear, as at 
large detuning the sideband essentially vanishes. Moreover, 
the strong asymmetry in xqd {(jj) at low temperatures has for 
consequence that the persistence of the peak is small for neg- 
ative detuning (see panel (a)). Finally, in Fig. [3l we show the 
influence of the temperature on the emission spectrum S{uj). 
The phonon sidebands grow with temperature and result in 
an increased emission through the cavity mode. 

This result suggests that the enhanced emission at the cav- 
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uj — u>o{meV) 





FIG. 2: Plot of the emission spectrum of the QD-cavity sys- 
tem in the presence of phonons (solid line) for different de- 
tunings D — ujo — oJc- a) D — —ImeV, h) D = OmeV, c) 
D = ImeV, d) D = 2meV (fcflr=10K). Comparison with no 
phonons (dashed line) is also given. 



ity mode, observed in several recent experiments, could be 
attributed to the phonon sideband mechanism. We point out 
however, that this cavity feeding phenomenon has been ob- 
served also when the exciton-cavity detuning is much larger 
than the typical width of the phonon broadbands, ie a few 
meV. These observations are accompanied by a superlinear 
dependence of the cavity mode PL on the excitation power. 
The phonon sideband mechanism, on the other hand, is ex- 
pected to provide a spectrum that depends linearly on the 
excitation power. The phonon sideband model is thus ex- 
pected to hold mostly at small detuning. The observations of 
cavity feeding at larger detuning are most likely due to multi- 
exciton emission, partially involving wetting layer states, as 
has been recently discussed [20l. [2TI. 143. [4^ . 
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uj — LOoimeV) 



FIG. 3: Emission spectrum of the QD-cavity system plotted 
for different values oiksT (arbitrary units), uic — ^oo = 1 meV . 



C. Influence of neighboring QDs 



the strong coupling by compensating for the cavity losses - 
a mechanism that can be traced back to the gain produced 
by the excitation of the additional QDs in the cavity. This 
enhancement cannot be reproduced by our model, as it would 
require accounting for nonlinear optical response. Eq. (|28p 
thus can model the presence of spectator QDs only in the 
limit of very small average population of each QD, for which 
the linear assumption holds. 



One major assumption of this model is that there is only 
one QD dot located within the region in which the cavity 
mode is extending. Given the density of the QD ensemble 
and the spatial extension of the nanocavity mode, it might 
well be that spectator QDs - i.e. additional QDs present in 
the cavity - contribute to the emission spectrum. These QDs 
are most likely weakly coupled to the cavity mode because of 
strong energy detuning or of smaller spatial overlap with the 
mode wave function. It has been recently suggested [isl . 
that, if these QDs are excited in addition to the main QD, the 
resulting emission spectrum is substantially modified, some- 
times even leading to a recovery of strong coupling in a sit- 
uation that would be of weak coupling if only the main QD 
was excited. In Refs. and [4^ . this effect has been mod- 
eled by an additional pump term acting on the cavity mode. 
Here, we can account for the presence of additional QDs in a 
natural way, by generalizing the expression for the single-QD 
susceptibility (|9]). The new susceptibility then reads 



XQD Ir,r ,tj) 



2 / 1 
^010 
" I 



where 



* (r) ** (r') XQD {io) + ^ *^ (r) ** (r') Xj 



(?8) 



(29) 



Here, the j-th QD has parameters defined analogously to 
those of the main QD. Starting from this expression, the 
derivation of the emission spectrum can be carried out anal- 
ogously to the single-QD case. In particular, when deter- 
mining the input field Qo {r, uj) as described in Appendix A, 
the virtual oscillating dipole will consist of a sum of terms 
originating from the different QDs, with relative weights Bj 
that express the contribution of each QD to the initial state 
of the emission process. It should be pointed out however, 
that the present model is based on the linear response to the 
Maxwell field. In Ref. [4^ the additional pump term enhances 



CONCLUSION 



In conclusion, we have shown that the Greens function for- 
malism is a powerful tool to relate quantitatively the usual 
atom-CQED parameters to the description of any QD-cavity 
system. We also extended this formalism to a QD weakly 
coupled to LA phonons. Thus we underlined that the differ- 
ence of a QD to a simple two level system is greatly enhanced 
when the quantum dot is placed inside a nanocavity. It makes 
possible that a PL peak of considerable amplitude remains at 
cavity frequency even for large detuning compared to the Rabi 
splitting, but limited to a few mev. This shows that the broad 
spectral features provided by the environment of a QD play 
a key role in the cavity-QD systems response. 
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Appendix A: Determination of the bare cavity 
electric field 



To determine Qo{r,ijj), we will follow the approach pro- 
posed in Ref. [2^. We have 



with 



V A V A (r, a;) - -Y^{r)£o {r, uj) = 



(Al) 



That is, with Ae(r) = e(r) - es, 



2 2 

V A V A {r, oj) - '^egSo (r, a;) = '^Ae{r)£o (r, uj) 



We define £s(r, oj) as a solution of 



V AV A£B{r,io)- -^esSB {r, uj) = , 



which is indeed a plane wave. At this point, we will use the 
using the virtual oscillating dipole method [3j|. We replace 
£b by a point source centered in ro: £b (f, i^) = BS{r — ro). 
Using the background Green's function defined as: 



V A V A ^13 {r, oj) - -^eBQB [r, oj) = S{r - r') , 



we have: 



£o (r, oj) = £b {r,ijj) 

dr'gB{r,r' ,Lo) — Ae{r')£B 



B 



S{r - ro) + GBir, ro, a;) — Ae(ro) 



(A2) 
(A3) 

(A4) 



In this expression, C/s(r, r',cj) is a slowly varying and espe- 
cially non-resonant function of lo. Then so does £o (r, oj) and 
finally: 

Qoir,Lu)^ y^£oir,Lj), (A5) 
with Qo {r, io) a function with no resonance. 



